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ABSTRACT 



This study was conducted to examine the eflects of discrete precipitation events on 
the short-term and seasonal evolution of ocean mixed la 3 'er temperature and salinity 
structure. This study was located at Ocean Station "P" (50°N, 145°W) in the Northeast 
Pacific Ocean. Two numerical modeling experiments were performed. The first was to 
simulate the response of the mixed layer to hypothetical discrete (isolated) precipitation 
events. This experiment showed that the effect of a single discrete rain event can var>' 
with season, with the mixed layer depth (MLD) at onset of the rain event, and with the 
strength of wind stress forcing. A single rain event can have lasting effects on mixed 
layer depth and temperature for up to 55 days after the event, depending upon the sea- 
son. 

The second experiment simulated quasi-realistic "complex" precipitation forcing, 
with a realistic distribution of synoptic events over a 13-month period. For this exper- 
iment, four different precipitation "intensities" w’ere tested. The values of temperature, 
salinity, and MLD predicted by the model were compared with observed Conductivity 
Temperature Depth measurements and with the values predicted using constant precip- 
itation forcing. In all experiments, the modelled MLD's approximated the observed 
MLD and temperature cycle. MLD's for all cases were too deep while temperature w'as 
estimated well in winter and w'as too cool in summer. Predicted salinity provided the 
greatest discrepancy between the modeled and the observed cycle. The 96 cm, 400-day 
quasi-realistic precipitation forced case best approximated the salinity observations 
though predicted salinity was fresher than observed in winter and saltier than observed 
in summer. Model results suggest that the amount of precipitation observed at Ocean 
Station "P" was too low’ to explain the observed and simulated ocean salinity and tem- 
perature structure for the year. 
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I. INTRODUCTION 



Salinity in the oceanic mixed layer and its change in depth with time can have a 
significant impact on upper ocean dynamics. Alterations of upper ocean salinity caused 
by precipitation and evaporation through buoyancy fluxes may also change the tem- 
perature profile of the ocean, altering the sound velocity. The purpose of this study w'as 
to examine the effects of discrete precipitation events on the evolution of the mixed layer 
depth (MLD) at Ocean Station "P" (50“N, 145“W). 

Previous one-dimensional mixed layer prediction experiments had failed to accu- 
rately explain and predict the behavior of salinity with time and depth. Possible reasons 
for this were; 

1. the amount of precipitation and evaporation was incorrectly specified, 

2. the precipitation and evaporation distribution with time, relative to wind events 
was incorrectly specified, and or 

3. the advection was incorrectly specified or neglected. 

This study focused specifically on the variation of precipitation and evaporation with 
time. The total amounts of precipitation and evaporation are also varied. Possible ef- 
fects of advection are discussed. 

A. LITERATURE REVIEW 

Kraus and Turner (1967) provided the fundamental turbulent kinetic energy (TKE) 
theoiy' on which many subsequent oceanic mixed layer models are based. This model 
considered the turbulent kinetic energy equation in a one-dimensional mixed layer 
scheme, using the approximately decoupled state of the mechanical and thermal energy 
equations. It was assumed that the heat input at the air-sea interface and the mass 
entrainment at the base of the mixed layer were uniform throughout the mixed layer for 
short time scales. Wind stress effects on the air-sea interface acted to increase the po- 
tential energy in the mixed layer by generating turbulent kinetic energy for mixing and 
entrainment. Entrainment was considered to be proportional to the wind energy input 
to the water column minus the work performed in overcoming the buoyant forces 
throughout the mixed layer. This type of turbulent erosion modeling considered 
buoyancy and wind forcing only, and it neglected viscous dissipation. This neglect of 
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dissipation has generally been considered a serious flaw. Kraus and Turner (1967) did 
not consider salinity effects in their mixed layer model. 

Subsequent works, such as Geisler and Kraus (1969), .Miropofskiy (1970), and 
Denman (1973) have attempted to consider dissipation by assuming dissipation to be a 
fixed portion of wind stress. The Denman (1973) model was used in the first mixed layer 
simulation experiment to use Ocean Station "P" data (Denman and Miyake, 1973). 

Miller (1976) was the first to include salinity and its effects on the density structure 
in a mixed layer model. Miller modified Denman (1973) to include salinity and used data 
from BOMEX Period III upper ocean soundings for periods of as long as 10 days for 
forcing and initial conditions. Four specific cases were examined by selecting distinctly 
different initial conditions for the salinity and temperature profiles and forcing data from 
BO.MEX surface observations. In three cases, .Miller examined the effects of three dif- 
ferent initial conditions while the fourth case included a precipitation event. 

The Miller study indicated that mixed layer salinity greater than the salinity of the 
water below the mixed layer may cause convective overturning due to density instability 
unless a sufficiently large temperature decrease occurred simultaneously across the 
interface. A large density jump across the interface acted to slow the mixed layer 
deepening rate due to an increase in potential energy. Without entrainment salinity flux 
effects, the cooling induced by entrainment was greater than the heating from the sur- 
face, resulting in net mixed layer cooling. With salinity included in the model, the layer 
deepened more slowly and entrained less cold water at the bottom, allowing surface 
heating to dominate and, therefore, the mixed layer temperature to increase. Overall, the 
results of these cases demonstrated the importance of salinity fluxes and the need to in- 
clude the salinity profile in the prediction of mixed layer deepening and resultant tem- 
perature structure. 

Miller demonstrated the potential significance of precipitation episodes in mixed 
layer modeling in a fourth case. His study demonstrated that precipitation can contrib- 
ute to salinity fluxes within the mixed layer. However, the fact that salinity regularly 
fluctuates by ±0.15% (due to horizontal advection) makes decreases in salinity due to 
rainfall difficult to detect in observations. Though Miller admittedly did not understand 
the nature of the effects of rainfall on the mixed layer, he did propose that even minor 
changes in the salinity flux due to precipitation may have a long term effect on the mixed 
layer thermal structure. 
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Paulus (1978) examined the effects of a salinity profile on density structure using a 
modified version of the Camp (1976) one-dimensional mixed layer model. The effects 
of a surface salinity flux on the mixed layer were examined using a constant evaporation 
and precipitation rate. Model initialization was accomplished with hydrocast data or 
XBT data with historical salinity profiles adjusted to realistically correspond to the 
temperature profiles. The incorporation of observed salinity profiles into the model had 
little effect on sea surface temperature (SST) and MLD during the summer regime (when 
MLD was shallow). However, when using the same forcing conditions for the winter 
regime the inclusion of a salinity profile acted to produce a warmer, shallower mixed 
layer than obtained when density was only a function of temperature. 

Paulus suggested that the retreat of the isothermal layer could be enhanced due to 
precipitation by slowing the deepening rate and shoaling tendency of the mixed layer 
under downward heat flux. With respect to salinity flux. Paulus concluded that imbal- 
ances in evaporation and precipitation had more pronounced effects on SST and .MLD 
in the long term than in the short term. His study concluded that this was a cumulative 
effect. Therefore, Paulus concluded that the addition of salinity effects to the mixed 
layer model was desirable. 

Garwood (1977) introduced a one-dimensional bulk model of the oceanic mixed 
layer that included dynamic instability with turbulent erosion. This model modulated 
the mean entrainment rate by diurnal heating and limited maximum layer depth by en- 
hancing the effect of viscous dissipation. This allowed a cyclic steady state to occur over 
an annual time period. It recognized the important effects of long term salinity fluxes 
on the mixed layer thermal profile due to surface heat fluxes, precipitation, and 
entrainment at the base of the mixed layer. This was the model of choice for this study. 

Garzon (1987) used an oceanic mixed layer model to determine the specific eflects 
of fresh water flux on the seasonal pycnocline in the mixed layer at Ocean Station "P." 
First, sensitivity of the model proposed by Garwood (1977) was examined by performing 
a hypothetical experiment using constant forcing with a wide range of precipitation and 
evaporation values. These results were compared to an identical model run with P-E 
(precipitation minus evaporation) = 0.0. This simple hypothetical case examined the 
relative importance of excess evaporation or precipitation over a time period of one- 
year. For the case with excess precipitation, a reduction in downward buoyancy and a 
decrease in surface salinity occurred, resulting in a shallower MLD than would be ex- 
pected if P-E = 0. For the case with excess evaporation, an increase in downward 
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buoyancy flux and an increase in surface salinity occurred. This yielded a deeper MLD 
than would occur with P-E = 0. 

For the annual simulation at Ocean Station "P," Garzon, using data from 1967, 
modeled temperature and salinity profiles with P-E = 0.0m 'year, P-E = 0.5m/year, and for 
monthly representative values of precipitation. Overall, MLD and surface salinity from 
the model well approximated those of the actual data. However, near the end of the 
study period, MLD and salinity values demonstrated increasing error particularly in ap- 
proximating significant anomalous "spikes." Garzon suggested that a more realistic 
representation of discrete precipitation and evaporation events would result in more ac- 
curate model output, especially for the long term. Additionally, Garzon adjusted model 
output for advection of salinity by subtracting the linear fit from the predicted salinity 
time series. 

Precipitation and evaporation events have been demonstrated to have a significant 
impact on mixed layer depth on both short and long time scales. These factors can act 
alone or may be coupled with such events as synoptic scale storms to modify the effects 
of these events. On a short time scale, daytime heating tended to cause shallowing of 
the oceanic mixed layer due to heating of the mixed layer despite the occurrence of 
evaporation which tended to promote deepening. Upward heat flux broke down the 
mixed layer causing deepening. Precipitation during this time period promoted shal- 
lowing. thus opposing this effect. These diurnal changes in mixed layer depth were 
moderated by annual effects that were associated with the annual radiation cycle. 
Strong winds, often associated with large scale synoptic events, when combined with low 
insolation values in winter, yielded deep mixed layer depths. However, precipitation 
events often accompanying these winds acted to lessen the effects of winds and were 
associated with discrete shallowing episodes that last for short periods of time. These 
efiects were cumulative. Just as a reduction or increase in wind speed in winter yielded 
a substantially shallower or deeper mixed layer depth in summer, a reduction or increase 
of rain during the same period yielded a deeper or shallower mixed layer depth in sum- 
mer. 



B. RESEARCH PLAN: LOCATION AND PERIOD OF STUDY 

For this study, a time period was chosen to provide the best overlap of available 
data bases. The starting date of 23 November 1968 with a study period of 400-days was 
selected for this study, since it provided data base overlap, that is, surface forcing and 



4 



CTD data was available for this time period. This numerical simulation of 400-days used 
a time step of one hour. 

The study area was between 49° X and 51° N, and 144° W and 146° \V centered at 
Ocean Station "P" (Fig. 1 on page 6). From 1949 to 1981 the waters around Ocean 
Station "P" have been occupied by ocean station vessels year-round. This provided nu- 
merous data bases, such as, XBT, CTD, and atmospheric observations for analysis. 

Tabata (1961,1965) described oceanographic and meteorological conditions at Sta- 
tion "P" in detail. The dominate pressure systems located in the North Pacific Ocean 
are the low pressure system or Aleutian Low in winter and the high pressure system or 
North Pacific High in summer. These two systems govern the distribution of mean 
surface winds over the Northeast Pacific Ocean. In winter the Aleutian Low directs air 
flow northward into the Gulf of Alaska, while in summer the North Pacific High directs 
winds southeast at Ocean Station "P." Additionally, monthly mean wind speeds are 
approximately twice as large in magnitude (13m's) in winter as in summer (6m s). 

Station "P" is located in a distinct oceanic regime (Fig. 2 on page 7). The surface 
waters of the region are under the influence of the weak westward flow. Specifically, this 
flow is composed of the West Wind Drift and the recirculated waters from the Alaska 
Stream. The former carries mixed waters of the Kuroshio and the Oyashio, which are 
generally indistinguishable by the time that they reach Ocean Station "P." To the east 
of this area, the current turns north, forming the Alaskan current. 

The vertical structure of this water is characterized by three zones: (1) the upper 
"seasonal" zone (O-lOOm), (2) the halocline (100-200m), and (3) the lower zone (> 200m). 
Within the upper zone the salinity is relatively low (~ 32.7ppm). With respect to sea- 
sonal effects, in winter the temperature and salinity in zone 1 are relatively homogeneous 
and isothermal. In summer a significant vertical salinity gradient exists in this upper- 
most zone, as well as a well defined thermocline, with a temperature decrease with depth 
of as much as 8°C per 20m. In the halocline, zone 2, salinity increases with depth by 
as much as l.Oppm per 100m, In zone 2 temperature decreases with depth typically, 
unless an inversion is present. In the lower zone salinity increases only gradually with 
depth, while temperature decreases gradually with depth (Fig. 3 on page 8). 

Tabata (1961) indicates that there is an annual average excess of precipitation over 
evaporation of ~ 0.5m. Garzon (1987) also estimates an excess annual precipitation. 
Despite this annual excess of precipitation, salinity within the water column does not 
decrease on an annual time scale, although it can vary as noted by Tabata (1961), 
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Fig. 1. Region Siinouiuling Ocean Station "P" (from Tabata, 1961) 




Fig. 2. Oceanic Regime at Station "P." (from Favorite, 1976) 
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Fig. 3. Example Temperature and Salinity Profiles 



8 



especially within the more dynamic mixed layer. Horizontal advection must account for 
this readjustment. As noted by Tabata, upwelling is weak and nearly insignificant within 
these waters, w'hen compared with horizontal advection. As a result this salinity balance 
appears to be maintained by horizontal advection as more saline water is transported 
into the area by ocean currents. 

Despite the consistent westerly current coming into this area, anomalous events 
have occurred resulting in horizontal transport of water into the area from a non- 
westerly direction. Tabata (1961) noted an unusual increase in temperature and decrease 
in salinity within the upper zone and halocline during spring 1957 through summer 1958. 
This was attributed to transport in w’ater lying south of Ocean Station "P" into the study 
area. The suggested cause for this transport was the shift of winds in the region associ- 
ated with unusual location of the atmospheric pressure system over the northeast Pacific 
Ocean. This displacement of the North Pacific High yielded westerly instead of the usual 
southeasterly winds, which could have generated southerly drift currents in the vicinity 
of the station. 

Overall, the currents in this area are weak. Except for occasional anomalous trans- 
port associated with unusual synoptic weather systems, salinity follows a roughly annual 
cycle maintained by horizontal advection. and the waters of the study area are relatively 
homogeneous horizontally. These factors combine to provide an excellent region for a 
one-dimensional modeling study of the thermodynamics of the upper ocean. 
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II. DATA SOURCES AND PREPARATION 



A. DATA SOURCES 

This study was performed using various sources. Atmospheric forcing data was or- 
iginally obtained from the National Weather Records Center, Asheville, North Carolina. 
This data included threc-hourly observations of wind direction and speed, air temper- 
ature, sea surface temperature, dew point temperature and cloud coverage in octaves. 
Atmospheric data was interpolated from three hourly to one hourly data for computa- 
tion of hourly surface fluxes. Precipitation forcing data was not directly available and 
was derived as explained later in this chapter. 

Two types of oceanic data were used. CTD/STD data was provided by the Marine 
Environmental Data Service, Ottawa, Canada. These salinity and temperature casts 
were digitized from depths of 0 to 1500 m at various intervals, with more concentrated 
data provided within the mixed layer. These casts were not taken at regular intervals, 
but were taken as operationally feasible. There were 119 casts available for analysis 
during the study period from station "P." Bathythermograph casts were not available 
for this time period. 

Three hourly atmospheric forcing data was used in this study to compute the fluxes 
of heat, radiation and momentum necessary' for forcing the mixed layer model with an 
hourly time step. Specifically, these boundar}' conditions were wind stress in two- 
dimensional components, incident solar radiation, heat loss (long and short wave), and 
precipitation and evaporation. The calculation of precipitation will be discussed later. 
Temperature and salinity profiles from STD taken 23 November 1969 were used to ini- 
tialize the model profiles, and CTD STD casts were used for verification. Bucket tem- 
peratures were used for comparison of with predicted temperature. 

B. VERTICAL INTERPOLATION 

The oceanic profiles, from both CTD/STD sources were vertically interpolated from 
0 to 199 m to a standard grid of 1 m spacing. Mixed layer depth, mixed layer temper- 
ature, and mixed layer salinity were derived from these interpolated profiles. There is 
not a standard way to compute mixed layer depth from a CTD/STD profile. For this 
study the mixed layer depth has been defined to be that depth at which the temperature 
profile changes exceed 0.10° C in 5 m. Mixed layer temperature and mixed layer salinity 
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values were obtained by averaging temperature from 1 m to the mixed layer depth as 
determined above. 
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III. THEORY 



A. OCEANIC MIXED LAYER STRUCTURE 

A familiarization with the definitions associated with the oceanic mixed layer ther- 
mal structure is required to understand the dynamics involved in mixed layer modeling. 
As defined in Garzon (1987) the oceanic mixed boundary' layer is assumed to be the fuUy 
turbulent region of the upper ocean bounded above by the air-sea interface and below 
by a relatively static water mass. The mixed layer is considered to be "quasi- 
homogeneous," that is, it is approximated to be well mixed. Thus temperature within 
the mixed layer is considered to be relatively constant. Fig. 4 on page 13 graphically 
depicts an oceanic vertical temperature profile. The turbulent mixed layer of depth h is 
contained within the upper most oceanic region, below which lies the turbulent 
entrainment zone, from depths z = -h to z = -h-d, with z positive upward. Within the 
entrainment zone there is a temperature jump of AT . Below the entrainment zone is the 
quasi-stable water mass and here temperature decreases with increasing depth. 

.Mixed layer structure (Fig. 5 on page 14) for salinity is similar to that for temper- 
ature (Fig. 4 on page 13). There is a relatively constant salinity within the turbulent 
mixed layer from z = 0 to z=-h, a salinity jump AS across the entrainment zone from 
z = -h to z = -h-d, and a gradually increasing salinity with depth below the mixed layer. 

An oceanic mi.xcd layer prediction scheme requires an accurate (or near accurate) 
assessment of the deepening or shallowing rate of that layer. Garwood (1977) accounted 
for both of these processes. The equations used by Garwood to best describe these 
processes are developed as follows. 

B. ENTRAINMENT 

The Garwood (1977) one-dimensional mixed layer model considers both turbulent 
erosion and dynamic instability. Fundamental to the solution of this model is an accu- 
rate assessment of the entrainment velocity, IV,. Starting with; 



4^ = IV - fV( - h) 

01 e V / 



(3.1) 



Assuming no mean vertical motion, W (-h) = 0, this equation becomes; 




(3.2) 
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Fig. 4. Example Temperature Profile 



The solution for IT, can be found by solving the total turbulent kinetic energy (TKE) 
budget written as: 



6i 2 



- V’W 



cU 

dz 



- VIV 



du 

dz 







-£(3.3) 



where: 

£ = is the Total TKE. 



(see Table 1 on page 20 for definitions of constants) 

The five terms in Equation (3.3) from left to right are: 
1. the time rate of change of turbulent kinetic energy, 
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Fig. 5. Example Salinity Profile 

2. shear production, 

3. buoyancy flux, 

4. turbulent dilFusion (or redistribution), and 

5. viscous dissipation. 

This equation may be examined for the entrainment zone. Within the entrainment zone, 
dissipation and shear production arc assumed negligible. This produces a mean turbu- 
lent kinetic energy equation for the entrainment zone of the form 
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across the entrainment zone of thickness <5. Integrating equation (3.4) yields: 



' ^ d F f~* 

-^■^dz=\ agTlV'dz- 

2 J . . 

-h-6 



f p' -i-h 

"2^-k 



T 

J-h-6 



(3.5) 



Assuming steady state, the time rate of change of turbulent energy may be neglected. 
Equation (3.5) becomes: 




Solving this integral produces equation (3.7). 



(3.6) 



^irn'’\_, + riv\ + I :!-< = o (3.7) 

When specifically examining this equation at the base of the mixed layer (or the top of 
the entrainment zone): 

‘rgTWL,^-^ir(f + -^)l_, (3.8) 



Kraus and Turner simplified the solution at this point by neglecting the redistribution 
term and combining the entrainment mean kinetic energy equation with the steady state 
turbulent kinetic energy equation for the mixed layer. Garwood does not neglect this 
term, but estimates its magnitude with the following parameterization: 



6z " ^ 2 + Po T, 



(3.9) 



Thus the redistribution term is equal to the mean turbulent energy available near the 
base of the mixed layer divided by the time scale required to transport this energy into 
the area. Parameterizing this time scale as 



T 



e 



_h 

\W'\ ’ 



(3.10) 



where | IP | is the root-mean square vertical turbulent velocity, I W | = (IP^)'-'^, and the 
time scale is approximated as the absolute value of the mixed layer depth divided by the 
absolute value of the vertical velocity. 
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Using Equations (3.8 to 3.10) and assuming that TW'\_^ = , 



g{yAT-p^S)\V,= 



< £ > I ir I 

h 



(3.11) 



Here vertical transport of turbulence exactly balances entrainment buoyancy flux. This 
equation may be solved for the entrainment velocity: 



IV \^y\<E> 

^ hg{(xAT-^AS) 



(3.12) 



To solve this equation for W, the values of < £ >, I W\ , T, h, and AT must be 
known while g and a are constants (Table 1 on page 20). From Garwood (1977) the 
total turbulent kinetic energy is solved by using: 



-§-{h<E>) = Go+G,-Bo-B,-D 

O I 



(3.13) 



In order to solve equation (3.13), the following equations for shear production at 
z = 0 and z=-h, and buoyancy damping 'production at z = 0 and z=-h, and dissipation 
are: 

wind shear production, 



<^o 




entrainment shear production. 



G/, = [ < z/ >^ + < V IV^, 



(3.14) 



(3.15) 



where the x and y components of the mixed layer current are: 



buoyant damping/produclion due to surface buoyancy, 




(3.16) 



buoyant damping due to entrainment, 



Bff = oighATlVg, and 



(3.17) 



dissipation. 



D = 2<E 



(3.18) 



C. MIXED LAYER THERMODYNAMICS: TEMPERATURE AND SALINITY 
BUDGETS 

Wind stress generated turbulence at the air-sea interface provides the energy to de- 
stabilize and entrain the water below. This mixing process deepens the mixed layer and 
mixes temperature and salinity between the surface and the entrainment zone. Thus, the 
conservation of heat requires that: 



Vertically integrating Equation (3.19) across the entrainment zone, z = — h — S to 
z = — h, and assuming V IV vanishes below z = — h — 5 gives 



AT =T{z = -h)- T{z = -h- S). 

Vertically integrating Equation (3.19) across the entire mixed layer, z=-h to z = 0 
gives 



Using Equation (3.20) and assuming T to be well mixed in z. Equation (3.21) becomes 



- AT Wg= T]V{-h) 



(3.20) 



where: 




(3.21) 
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er rff’io) TW{-h) 

Cl ~ h ^ h 

or 

ST _ Qo 

Cl pCph h 

Here : 

1. Af=T-T{-h-d), and 

2 . QoIpC, = -TW{0). 

The one-dimensional salinity budget for the oceanic mixed layer is: 

ds _ sirs’ 

8i cz 



(3.22) 



(3.23) 



(3.24) 



It is assumed that: 

1. S(z) is well mixed from the bottom of the mixed layer to the top of the mixed layer, 

2. S(z) has a near discontinuity between z = -h and z= -h-(), and 

3. the vertical flux of S (or S’ 11") vanishes at z < —h — b. 

Fig. 5 on page 14 depicts this condition. Here the turbulent region within the mixed 
layer yields a near homogeneous state with respect to salinity. Examining Equation 
(3.24) and integrating gives: 



(>-W) — 

8S 



dz = 



>-h 



-h-6 



clV'S' 
6z 



dz 



(3.25) 



With h variable in time, application of the Leibnitz Rule of integration and assuming 5 
is vanishing small. 



AS 4 ^ = - irS'i -h) 4- 1F'S'( -h-b). (3.26) 

c / 

Assuming that turbulence is negligible below the mixed layer, i.e. 1F'S'( — — <5) ~ 0. 
This makes the "jump condition" for salinity entrainment into the mixed layer: 

AS 4^ = - TF^( -h). (3.27) 

Cl 



IS 



where: 

AS = S{z = - h) -S{z - h - S) 



D. SALINITY EFFECT ON MIXED LAYER DYNAMICS 
1. The Shallowing Effect 

Kraus and Turner's (1967) model provides an approximation for the non- 
entraining or shallowing mixed layer, IF, = 0, 



C, 

^ o.gQol{pCp) 



= L. 



(3.28) 



Here L is the Obokov length scale and equation (3.28) is the solution when there is no 
entrainment. The surface buoyancy flux is = {o.gq^j{pC^) and U. is the surface friction 
velocity, with bl = {pJp)C^U'^ . The value of h is inversely proportional to 2o- The 
mixed layer temperature is even more sensitive to than is h; 

= Q^lipCph) ^ "T (3.29) 

Cl /i 



If the buoyancy flux attributable to precipitation and evaporation is included, Equation 
(3.28) becomes 

C, ul 

/? = -r . (3.30) 

^2 o.gQ,i{pCp) + pg{P-E)S 

Where the salinity concentration coeflicient, = 0J6x\0-^Kgjgm , and P-E is the net 
precipitation minus evaporation in m sec. With the salinity flux included, the surface 
buoyancy flux becomes; 

^0 = ^-SQaipCp + h{P - E)S, (3.3 1) 

and the effect on mixed layer temperature becomes: 

= (QolpC,)lhoiQo{Qo + mP - E)S)- (3-32) 

C I 



As a result, precipitation minus evaporation influences mixed layer temperature indi- 
rectly. For positive go (r^^t the surface), a positive (P-E) will decrease h, con- 

centrating the heat into a shallow layer. For the case where evaporation exceeds 
precipitation a decrease of mixed layer temperature with respect to time results. 
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2. The Deepening Effect 

For entraining, deepening mixed layers, the buoyancy discontinuity at the bot- 
tom of the mixed layer must be augmented to include a possible salinity jump: 

Ab = cxgAT - ligAS (3.33) 



where: 



AS = S -S{-h-6) , 



S is the mixed layer salinity, and 

S( — — 5) is the salinity just below the mixed layer. 

The mixed layer salinity budget is: 

(3.34, 

ct h h 

Equation (3.34) demonstrates that the change in salinity with respect to time 
depends on fresh water flux at the surface, entrainment velocity, salinity advection, 
mixed layer depth, and upon the salinity jump at the base of the mixed layer. 



Table 1. MODEL CONSTANTS 



Model Constants 


constant 


definition 


value 


a 


thermal expansion coefficient 


~ 10- A’-' 


p 


salinity expansion coefficient 


~ 8 X 


Q 


tuning coefficient 


~ 1 


Q 


tuning coefficient 


^ .2 


Cr 


specific heat 


~ 40007 . kg-' . K-' 


p 


sea water density 


~ 1026. 
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IV. HYPOTHETICAL EXPERIMENTS 



A. EFFECTS OF DISCRETE PRECIPITATION EVENTS 

The purpose of this study was to determine the lasting effects of discrete precipi- 
tation events on the mixed layer. For "case 1" in this study short-term precipitation 
events of 10 cm in one day with no evaporation was applied monthly to the Garwood 
(1977) model over a 365 day study period. This yielded twelve separate events for de- 
tailed study. At all other times during the year both evaporation and precipitation were 
assumed to zero. The resulting output from the model was examined individually for 
30 day periods following each rain event to better reveal long term trends. The precipi- 
tation event was applied on the fourth day of each monthly period. Actual temperature 
and salinity profiles taken from a CTD cast for November 23, 1968 were used as initial 
conditions. Heat and momentum fluxes were computed from the actual surface 
meteorological and observations. 

Fig. 6 on page 22 and Fig. 7 on page 23 are a conglomeration of the results of these 
twelve separate events. Other figures in this chapter represent the individual events that 
were selected for closer examination at a 60 day period. 

The effects of precipitation events on the M LD are graphically represented by the 
figures in this chapter. The equations developed in Chapter 3 will be used to explain the 
results and are repeated here: 



dT _ Qo 

dt pCph h 



(4.1) 



The change in mixed layer temperature (T) with respect to time is a function of: 

1. MLD, h, 

2. entrainment rate, IF, 

3. temperature jump at the base of the mixed layer, AT, and 

4. net downward heat flux, Qq. 

For a case of no entrainm.ent, IF, = 0, the temperature will rise for > 0, and it will fall 
for Qq < 0. The rate of rise and fall is inversely proportional to h. If ^ 0 and there 
is entrainment, IF, > 0, the mixed layer temperature can only decrease if AT > 0. The 
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IZ-P = 0, (c) winJ speed, and (I) net suil'aee heat flux. 



->■> 




Fig. 7. Coiiglonieralion of Six Events (Jun to Nov 1969); (a) didcrciuial MLD, 

(b) didcrciuial temperature, (c) did'erential salinity, (d) MLD witli E-P = 0, 

(c) wind speed, and (f) net surface heat Ilux. 
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rate of decrease is greatest for small h. Various factors effect U\, and will be considered 
later. 

Two model predictions may be compared. In this experiment the first case included 
the discrete precipitation event, while the second case had no events, with E-P = 0 for the 
entire model integration. When subtracting the results of one study using equation (4.1) 
from another such study for the same time using the same equation, the results provide: 

1. differential .MLD (/z, — /zj), 

2. differential temperature (7, — Tj), and 

3. differential salinity (S, — Sj). 

Differentiating the differential temperature with respect to time yields: 



Therefore, in order to provide a value other than zero for ^{T^ — T^jct , the differential 
temperature, one or more of the following conditions must exist: 

1 . /Z, /z,. 

2- el ^ fi f2i ^ntl 



Again, considering the many variables that may alter the entrainment rate, there are 
numerous factors which may effect this differential temperature, making interpretation 
of even a simple precipitation event complex, but more easily exantined than the case 
with a more realistic and complex E-P time series of forcing. 

Initially, in all twelve events with all forcing values equal (for the first five days), the 
differential temperature (57= 7, — Tj must equal zero. This can be seen in Fig. 6 on 
page 22 and Fig. 7 on page 23 for the first five days of each 30-day period. Upon initi- 
ation of the precipitation event at day five of each events h^ = and A7, = A7j , but 
IF, is instantaneously altered in the precipitation case. In this case with E-P < 0,there 
results in an immediate increase of surface buoyancy, , 




(4.2) 



3. A7,#A72. 



B^ = o.g -0gS{E-P). 



An increase in acts to decrease the total TKE from equation (3.13), 
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{h< E > ) — Gq + Gf^ — Bq — — D 



and, therefore, entrainment rate from equation (3.12), 



" gh{aAT-pAS) ' 

As IF, decreases, so will MLD, which acts dynamically to shallow the MLD relative 
to the case with P-E = 0. Thus a precipitation event increases the downward surface 
buoyancy through decreasing salinity, which ultimately results in a shallower MLD. 
This initial result should be apparent throughout the twelve experiments where there is 
precipitation. 

While the immediate change in MLD through reduced entrainment is expected on 
a short time scale (during the rain event), longer term changes in mixed layer dynamics 
may result from a single short-term precipitation event. In response to mixed layer 
shallowing, MLD and the temperature jump at the base of the mixed layer will change 
from that of the non-precipitation case. At this point all three possible conditions re- 
quired to yield an increase in the differential temperature may exist. 

Thus analysis of the effects of a discrete precipitation event requires a case by case 
examination. Any differences in MLD will yield a non-zero magnitude of I/i, — /jjl 
which will effect the differential temperature. A fluctuation in differential MLD will also 
moderate the effects of Qa , AT^ — ATj, and fF,, — If ',2 . Any difference in the temper- 
ature jumps, IAT, — ATjl 7^0 , will also jield an associated change in the differential 
temperature and a transient response, as well as the altered value of fF, . 

B. ANALYSIS OF PRECIPITATION EVENT EFFECTS 
I. Event #I (Dec 1968) 

Fig. 6 on page 22 represents the MLD output for the model with E-P = 0. This 
figure indicates that this event began in the late-fall season when the mixed layer was 
at 80 m and deepening. Fig. 6 on page 22 and Fig. 8 on page 26 indicates that there 
was a net heat loss {Qa < 0) during this period associated with a mid-latitude, winter re- 
gime. There were moderate winds during this time period {U = 5.0ml s) , and maximum 
values of = 20.0w/s which occurred immediately before, during and after initiation 
of the precipitation event of day 5. 
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Prior to initiation of the precipitation event all differential values equaled zero 
(Fig. 6 on page 22 and Fig. 7 on page 23). This was true for all twelve events. Upon 
onset of precipitation at day four, MLD decreased from 65m to 10m, yielding a decrease 
in the differential MLD. This was coincident with shallowing of the mixed layer asso- 
ciated with a lull in wind. Approximately 12 hours after precipitation onset, however a 
moderate wind event followed immediately by another of magnitude U ~ S.Omjs oc- 
curred. These increases in wind stress occurred during the precipitation event and re- 
sulted in two increases of the differential xMLD and temperature (points a and b). Thus, 
the layer deepened, due to decreased causing the layer to cool. 

Immediately following events a and b the wind decreased to near zero causing 
a large decrease in .MLD. Here the differential .MLD reached a maximum, as the .MLD 
in case 1 shallowed to near the surface. Associated with this decrease in .MLD in case 
1, was a negative differential temperature. This was the result of the shallower .MLD in 
case 1 allowing for the increased heat loss by the layer with ^o>0 (net heat loss) in 
winter (point c). The differential .MLD increased then decreased, as induced by the 
onset and drop-off of a wind event (points d and e). Again this produced an associated 
positive and negative differential temperature, respectively. At day six a strong storm 
occurred, (L’ ~ I3.0mjs) (event f) followed immediately by a series of larger fluctuations 
in wind speed. These storms brought the MLD of case 1 to near that of case 2. In- 
creased buoyancy, causing decreased entrainment, caused the MLD in case 1 to deepen 
less than that of case 2 accounting for this significant differential temperature event at 
day six as case 1. This entrained less cold water than in case 2. After day six a series 
of winds (U ~ 5.0ml s) continued to mix the layer, bringing the differential .MLD's and 
temperatures closer to zero. Although the effect of the precipitation event on the mixed 
layer decreased with time, the layer temperature and depth were effected measureably for 
this experiment period of 60 days. 

There was an immediate decrease in salinity as precipitation diluted the layer, 
and this decrease continued until cessation of the precipitation event. Mixing by 
entrainment decreased the differential salinity to a value of -0.4ppm, after which time it 
rapidly approached that of case 2, aided by wind mixing. 

2. Event #2 (Jan 1969) 

As in Dec 1968, prior to the precipitation event the MLD was relatively deep 
{h ~ 90m) and continued to deepen during this winter regime. The values for still 
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indicated a net loss to the atmosphere and strong winds were predominate 
{U ~ 5.0ml m) . 

After the onset of the precipitation event, the results of the event were ver>^ 
similar to those of event ^1. Therefore, no specific figures were provided other than 
those found between days 30 and 60 in Fig. 6 on page 22. 

At the beginning of the precipitation, differential .MLD generated a large de- 
crease ( ~ — 100.0/m) as the MLD of case 1 nearly reached the ocean surface. This was 
concurrent with a wind speed lull for two days, enhancing this .MLD shallowing. Asso- 
ciated with this .MLD decrease was a decrease in the differential temperature as the 
shallower .MLD in case 1 enhanced the temperature decrease associated with a net loss 
of heat in winter, < 0. By day 40 a series of strong wind events ((/ ~ 1.5ml s) acted 
to converge the .MLD of case 1 back to near that of case 2. Likewise, the differential 
temperature approached zero in a stepwise fashion as the temperature in case 1 ap- 
proached that of case 2. Decreased entrainment rate due to increased buoyancy did not 
allow the MLD of case 1 to remain as deep as that of case 2. This caused the temper- 
ature of case 1 at day 23 to exceed that of case 2. It did not, however, exceed the in- 
crease seen in Event #1 as the difference in the MLD's between case 1 and case 2 were 
smaller in this event. 

There was an immediate decrease in salinity as precipitation diluted the layer 
and this decrease continued until the end of the rain event. .Mixing by diffusion and 
advection decreased the differential salinity to a value of -0.3ppm, after which time it 
rapidly approached that of case 2. 

3. Event #3 (Feb 1969) 

The .MLD prior to the precipitation event was deep (h ~ 95//i) and deepening 
during the winter regime. The value of still indicated a net loss of heat to the at- 
mosphere, but the magnitudes of Qg and of the wind speeds were smaller than for the 
previous two events. 

After onset of the precipitation event in case 1, the results of the event were very 
similar to those of Event #1 and #2. As in Event #2, no specific figures were provided 
other than those found between days 60 and 90 in Fig. 6 on page 22. As the rain began, 
differential MLD generated a large negative event (~ — 70.0//i) as the .MLD in case 1 
nearly shallowed to the ocean surface. This was concurrent with a wind lull of four days, 
enhancing this shallowing. Associated with this MLD decrease was a decrease in the 
differential temperature as the shallower .MLD in case 1 increased the temperature drop 
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associated with a net surface heat loss. During the end of the precipitation event there 
was a small peak in both the differential MLD's and temperatures associated with a wind 
speed maximum. By day 70 a series of strong wind events {U 5.0ml s) caused the .MLD 

of case 1 to converge to near that of case 2. Likewise, the differential temperature ap- 
proached zero in a stepwise fashion as the temperature in case 1 approached that of case 
2. Decreased entrainment rate due to increased buoyancy did not allow the MLD of 
case 1 to remain as deep as that of case 2. This caused temperature of case 1 at day 70 
to exceed that of case 2. It did not, however, exceed the increase seen in Events #1 and 
#2, as the difference in the MLD's between case 1 and case 2 were smaller in this event. 

There was an immediate decrease in salinity in the first case, as in the first two 
events, as the precipitation diluted the layer. This decrease continued until cessation 
of the precipitation event. Mixing by entrainment decreased the differential salinity to 
a value of -0.3ppm, after which it rapidly approached that of case 2. 

4. Event #4 (Mar 1969) 

Examining the MLD model output for the E-P = 0 case (Fig. 6 on page 22), it 
appears that this was the period of deepest MLD (/i > 100m). This was also the begin- 
ning of the period of transition from the colder winter regime to the warmer summer 
regime. Here the MLD began to shallow due to decreasing wind speed and Qq > 0. 

This event provided different results than Events #1 to #3, and Fig. 9 on page 
30 is presented for better resolution. As in the previous events, there was an initial de- 
crease in both the differential MLD and in temperature (point a). This was associated 
with the precipitation event and a large decrease in U from 18.0 to near 0.0 m's. During 
this shallowing period there was a slight wind increase causing U to peak to 5.0 m/s, and 
a deepening in both case 1 and case 2. This was less in case 1 than in case 2. This dif- 
ferential in MLD generated a positive differential temperature (point b). The positive 
differential temperature immediately decreased as the wind dropped to zero, and the 
MLD became very small, order 1 m (point c). Net heat loss cooled the upper layer in 
case 1 more than in case 2, generating a severe negative spike in the differential tem- 
perature (~ — 0.10°C). A series of mild wind speeds until day 98 yielded various negative 
and positive peaks in the differential MLD's and temperatures. 

At day 100 a strong wind (1/ ~ 8.0m/s) returned the MLD in case 1 equal to that 
of case 2, making the differential MLD zero (point d). This mixing decreased the dif- 
ferential temperature to near zero, but it remained slightly colder as the entrained water 
from below the mixed layer remained warmer than water of the previous shallow layer. 
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Strong winds at days 113, 130, and 140 continued entrainment, thus, diffusing the colder 
layer in case 1 with the warmer water below, and bringing the differential temperature 
closer to zero (points e to g). 

Near the end o^this 60-day experiment there was a very strong mixing event at 
day 147, induced by a U peak of 15.0 m/s with a duration of 3 days. Almost 60 days 
after the single precipitation event, there was still an effect on the mixed layer. Exam- 
ining the differential MLD and temperature, it appeared that the MLD in case 1 did not 
deepen as much as in case 2. This deeper MLD caused the temperature in case 1 to begin 
to exceed that of case 2 (point h). 

Following the cessation of precipitation, in the first case, diluted salinity was 
mixed by entrainment, decreasing the differential sahnity to a value of -0.4ppm. After 
this it rapidly approached that of case 2 by wind mixing. 

5. Event #5 (Apr 1969) 

Near the end of the March study period (Event ^4) the winds decreased and the 
MLD shallowed (Fig. 6 on page 22). After the wind speed lull from day 123 to day 127, 
there was a steady rise in wind to the peak value at day 147 of L’ ~ lA.Omjs . This 
gradual increase provided the TKE to mix the layer to to the deepest point for the year 
{h ~ 125.0m). This was followed by a period of shallowing, induced by an abrupt fall in 
wind speed (Fig. 10 on page 32). The value of positive and increasing, net heat 

flux. This served to provided a net heat gain in this summer regime, as opposed to the 
heat loss that occurred with the winter regime. 

Fig. 10 on page 32 demonstrated the initial decrease of MLD in case 1, to be 
greater than that of case 2 (point a). Due to the low wind speeds coincident with the 
initiation of the precipitation event, the mixed layer of case 2 was also relatively shallow 
{h ~ 21.0m) . Thus the increase in the differential temperature was less than that of one 
day later (at point b) following a wind-driven deepening event. This produced a large 
decrease in the differential MLD. A series of shallowing and deepening events followed 
and were associated with differential temperature adjustments. This brought the differ- 
ential MLD and the differential temperature back toward zero. 

The largest decrease in differential salinity in Fig. 6 on page 22 occurred in this 
event. In this event the precipitation event was initiated after a long period of low 
winds, providing a shallow MLD prior to onset of the precipitation event. This may 
have been an unrealistic example. Given that increased wind speeds generally ac- 
company large precipitation events like that of day 124, this hypothetical event would 
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have been unlikely to occur, but it permits an interesting examination of model dynam- 
ics. 

6. Event #6 (May 1969) 

Between day 150 and day 180 the MLD in case 2 was similar to that of Event 
#5, but more shallow. This shallowing was the result of decreased wind forcing as seen 
in Fig. 6 on page 22. Net heat gain to the ocean surface continued to increase as sum- 
mer approached, and some of the largest net downward heat fluxes occurred during this 
period (<2omax ~ ^OOflO- 

The dynamics of this event resembled those of Event #5. No 60-day period fig- 
ures were provided. Fig. 6 on page 22 provides adequate resolution. As before, there 
was an initial shallowing producing a positive differential temperature associated with 
the shallower MLD and the summer heating regime. Until day 165, the differential 
MLD and the differential temperature resembled those of the previous events. After day 
165, however, the resulting differential temperature profile differed from those of the 
previous events. Here, the differential temperature reached a post-precipitation mini- 
mum value of 0.05°C and continued to increase slowly through the rest of the study 
period. 

This steady increase in differential temperature, when compared with the tem- 
perature plot of the mixed layer, indictated that the mixed layer was warming faster for 
case 1 than case 2. The .MLD's were equal in depth during this period, so that this 
condition could not have been due to case 1 being shallower than case 2. The possible 
causes were decreased IF,, due to decreased salinity, and/or increased AT. An increased 
AT could be caused by the mixed layer in case 1 being warmer than that of case 2, and 
the temperature of the water mass below being the same for both cases. If so, in the 
summer regime this steady increase in the difl'erential temperature would appear to have 
been "self-perpetuating." 

There was an immediate decrease in salinity in the first case, as in the other 
events, as the precipitation diluted the layer. This decrease continued until cessation 
of the precipitation event at day 156. Mixing by entrainment decreased the differential 
salinity to a value of -0.3ppm, after which time it rapidly approached that of case 2. 

7. Event #7 (Jun 1969) 

Though the MLD for case 2 began relatively deep {h ~ 70m) it decreased to the 
shallowest mean period value {h ~ 20m). This was caused by an initial wind event at day 
181 {U ~6mjs) , and followed by a relative lull from day 182 to 205. Another wind 
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increase occurred, similar in magnitude to the first at day 205 (Fig. 1 1 on page 34). The 
net heat flux remained strongly downward, but not as strong as in the previous study 
(Fig. 7 on page 23). 

Initially, the MLD for case 2 was shallow (/ij ~ llm) . At the beginning of 
precipitation, the MLD of case 1 shallowed to only 3m, yet there remained a small dif- 
ferential MLD, producing a small differential temperature <0.1°C (Fig. 11 on page 
34,point a). This difference decreased with continuing low level winds which acted to 
deepen the MLD in both cases and yielded a decrease in the differential temperature to 
near zero (point b). 

The most significant deepening event of this period then occurred (point c), 
which produced the largest differential temperature of the period ( ~ 0.5“C) . The mixed 
layer of the rain case then behaved as expected, stepping differential MLD and differ- 
ential temperature back toward zero. The differential temperature for this period was 
much larger than for the other periods studied. This appeared to be due to the very' 
shallow MLD's achieved by case 1, enhancing the effect of the downward summer heat 
flux. 

The shallow MLD at the onset of precipitation did yield a stronger than normal 
differential salinity ( ~ — O.OlOppm). The differential salinity acted as in the previous 
events, returning to an equilibrium value of ~ — O.Olppm . 

8. Event #8 (Jul 1969) 

This event coincided with a shallow (/i ~ 30//i) trend in the MLD for case 2. 
This was associated with low wind speeds {U < 5m/s) (Fig. 7 on page 23 and Fig. 12 on 
page 36). Xet heat flux appeared to be steady ~ 180 IF) . 

There was an initial positive increase in differential temperature which occurred 
in a stepwise manner, commensurate with the differential MLD (Fig. 12 on page 36). 
This occurred until wind speed decreased to zero at day 216 (point a). From day 216 
to 226 (points a to b) the differential temperature remained nearly constant at 
~ 0.065® C, with only a slight increase in time. The .MLD in case 1 equaled that of case 
2 (and the differential MLD ~ 0) in this period. This condition w’as similar to the slight 
increase which was noted in for event #6, and the physical explanation is the same. 

After day 226 (point b) there was a large increase in the differential MLD and 
differential temperature due to a U spike of ~ 4.2fnjs. This acted to cool the mixed layer 
of case 2 more than that of case 1 and yielded a positive increase in the differential 
temperature ( ~ 0.1®C). The differential temperature then fluctuated around 0.22°C for 
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Fig. 12. Event //8 (Jul 1968): (a) dillcrcntial MLD, (b) dinbrcntial temperature, 

(c) dill'erential salinity, and (d) wind speed 
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the remainder of this 60-day period as mild winds allowed shallowing of the mixed layers 
in cases 1 and 2 with intermittent weak spikes in wind speed producing some oscillation 
around this mean. 

There was an immediate decrease in salinity in the first case, as in the other 
events, as the precipitation diluted the layer and this decrease continued until cessation 
of the precipitation event. Mixing by entrainment decreased the differential salinity to 
a value of -0.3ppm, after which time it rapidly approached that of case 2. 

9. Event #9 (Aug 1969) 

With the exception of the initial small increase in differential temperature, this 
event produced results very' much like that of Event #8. This small fluctuation was ac- 
counted for by the weak wind mixing {U ~ Anils) occurring upon onset of precipitation. 

10. Event #10 (Sep 1969) 

During this study period the MED in case 2 was relatively shallow, and begin- 
ning to deepen, as there was a seasonal shift from summer to winter (Fig. 7 on page 
23). This was associated with steadily increasing wind speeds (Fig. 13 on page 38). As 
the seasonal regime began to shift to a winter regime, the shifted from one of net 
downward heating toward one of cooling, although, during this period net heating re- 
mained the slightly positive. 

As with the results of Event #9, these experimental results resembled those of 
Event #8 without the early fluctuating. In this event, prior to precipitation onset, the 
MED in both cases had shallowed to near surface due to a wind speed lull from day 273 
to 275 (Fig. 13 on page 38). A slight fluctuation at day 274 (point a) generated a small 
downward event in differential MED (~ — Im) which promoted a slight upward differ- 
ential temperature. With a follow on U event at day 275 (point b) of ~ Amis, the largest 
of all positive differential temperature in the twelve events was generated. This occurred 
along with a differential MED of ~ — 25m. Wind events followed which increasingly 
deepened the MED of case 1 and brought the differential temperature back to near zero, 
and constantly decreasing. 

This event produced the largest differential salinity of all, concurrent with this 
anomalous differential temperature fluctuation. Again, the veiy' shallow MED estab- 
lished before precipitation onset appeared to have contributed to this effect, and equi- 
librium was quickly reestablished at differential salinity ~ — O.Olppm. 
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11. Event #11 (Oct 1969) 

In this period the MLD of case 2 was becoming deep again (h ~ 60m) and 
continued to deepen in this early winter regime (Fig. 7 on page 23). Increasing wind 
speeds and near zero net heat values accounted for this effect. 

This return to the winter regime bidded a return to the differential MLD and 
differential temperature profiles seen in the first three events. Initial cooling of the 
shallower MLD of case 1 generated a negative spike in the differential temperature 
which moved positive as the decreased entrainment rate of case 2 inhibited entrainment 
of the layer. Further mixing acted to equalize MLD's and return the differential .MLD 
and differential temperature toward zero. After the differential temperature flux moved 
to a positive value it remained at a larger level than observed in the first three events, 
as it did in the summer regimes. This indicated that the mixed layer was in transition 
from the summer to wi ner regimes and MLD's were shallower than those of the first 
cases. 

There was an immediate decrease in salinity in the first case, as in the other ex- 
periments, as the precipitation diluted the layer and this decrease continued until ces- 
sation of the precipitation event. .Mixing by entrainment decreased the differential 
salinity to a value of -0.2ppm, after which time it rapidly approached that of case 2. 

12. Event #12 (Nov 1969) 

This event was very similar in forcing conditions and resulting output to that 
of Event #11. In this event the .MLD in case 2 was deeper than that of Event #11 as 
early winter deepening continued, and mean wind speeds appeared to be stronger with 
larger peak wind events (Fig. 7 on page 23). This deeper .MLD yielded a larger differ- 
ential .VI LD after precipitation onset, as well as a lesser positive peak of differential 
temperature. This also yielded a similar differential salinity perturbation. 

C. EXPERIMENT SUMMARY 

The results of these events indicate that the effects of discrete precipitation events 
varied greatly as a function of low frequency seasonal heat flux forcing and higher fre- 
quency wind forcing. These two forcing boundary conditions strongly influenced the 
mixed layer thermal characteristics, as well as the salinity. 

The average deeper, colder mixed layer in winter was significantly shallowed by a 
single precipitation event. This shallowing allowed for the increased loss of the thermal 
energy of the shallower mixed layer to the colder atmosphere, supporting a temperature 
decrease in the mixed layer. Wind events that followed precipitation deepened the mixed 
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layer, but it did not deepen as much, nor as quickly as in the case without precipitation. 
Follow-on wind events deepened the mixed layer until entrainment returned the MLD 
to that of the case without precipitation. Entrainment had a similar effect on the dif- 
ferential mixed layer temperature. 

The effects of the precipitation event in winter w^ere modulated by the discrete 
deepening and shallowing events that occurred just prior to precipitation onset. A 
deeper MLD increased the magnitudes of the differential values that occurred as a result 
of precipitation shallowing. A shallow^er MLD operated in the opposite sense to de- 
crease these values. Deepening and shallowing events occurring concurrent with pre- 
cipitation had a similar effect, acting to enhance or decreasing differential values. 

In summer the average shallow'er, warmer mixed layer was less significantly shal- 
lowed by these precipitation events. Given that the mixed layer was already relatively 
shallow, the decrease in MLD was not as significant as in winter. Significant thermal 
stratification in summer yielded a large increase in mixed layer temperature with precip- 
itation induced shallowing, as the thermal energy input at the top of the mixed layer was 
concentrated into a shallower mixed layer. This increased mixed layer temperature 
produced a more stable mixed layer as the temperature jump at the base of the mixed 
layer increases. The increase in temperature jump tended to further inhibit entrainment 
and slowed deepening of the mixed layer. Deepening slowed enough that the differential 
temperature actually increased with time. 

The discrete precipitation events could have had a long term effect on the mixed 
layer. This was especially true in summer when MLD's were shallower. In cases where 
strong winds generated strong mixing events, the mixed layer continued to demonstrate 
increased buo^'ancy, as long as 55 days after the precipitation event, as the differential 
MLD became negative and differential temperature yielded a positive fluctuation. 
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V. PROCEDURES, RESULTS, AND ANALYSIS 



A. PROCEDURES 

The 400-day mixed layer simulation, as described previously, was conducted to ex- 
amine the accumulative effect of a time series of rainfall events in explaining salinity and 
temperature changes with time and depth. Discrete daily precipitation data was not 
available for this study. Only monthly totals were available. As a result a method was 
devised to simulate daily precipitation values using the available oceanographic and 
meteorological observations. With the lack of surface pressure and upper atmospheric 
data, precipitation prediction became difficult, thus prediction of localized atmospheric 
convective activity was difiicult. Utilizing wind and cloud information, however, an at- 
tempt has been made to predict precipitation associated with frontal activity. This 
technique was based on meteorological principles developed by Petersen (1956). 

Specifically, this procedure utilized wind speed and direction, and cloud coverage 
data (Fig. 14 on page 42). With cloud cover data available in octaves, precipitation was 
allowed for cases with cloud coverage equal to or greater than 7/8. Considering pre- 
frontal and frontal rain, precipitation was allowed to occur for cases with wind directions 
from 000° (northerly winds) clockwise to 220° (southwesterly winds), and wind speeds 
greater than 10 kts. This accounted for the southerly to southwesterly winds that pro- 
ceed a front and the "confused" northerly to southeasterly winds that accompany frontal 
passage. For post-frontal rain, precipitation was allowed to occur for cases with wind 
direction 220° (southwesterly winds) clockwise to 000° (northerly winds), and for 
stronger wind speeds, greater than 25 kts. This accounted for the strong southwesterly 
to northwesterly winds that accompany a post-frontal regime and ensured consideration 
of all possible wind directions. Additionally, the post-frontal case required that the sea 
surface temperature be greater than the air temperature, as the cold frontal air passes 
over warmer water. Precipitation amount was adjusted to yield the 400-day precipitation 
total of the observed monthly data. The resulting rain with time is contained in Fig. 15 
on page 43. 

Tliis "complex" alogarithm for precipitation forcing scheme was applied to the 
OPEL model at a range of possible amounts that may represent normal interannual 
variability in the North Pacific. A "medium" precipitation level was chosen by adjusting 
total 400-day output to equal that of the summation of the 13 month period from 
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Fig. 14. Front and associated nind directions and isotherms 
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Fig. 15. Complex Rain Model Output and Winds for 400 day period: (a) rain 
model output for 400 days from 23 Nov 1968 adjusted to equal actual 
monthly totals, and (b) wind for same time period. 
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December 1968 to January' 1970 (Table 3 on page 45). The "low'" level was chosen to 
equal two-third's of the total 400-day precipitation total of the medium level case. The 
"high" level case was equal to one and one-half times the medium case and the "very 
high" case was made equal to three times the medium case (Fig. 16 on page 46). The 
same wind stress and heat flux atmospheric forcing values and initial conditions de- 
scribed in chapter II were used for all four model runs as in chapter IV. 



Table 2. PRECIPITATION TOTALS FOR THE 400-DAY STUDY PERIOD 



Precipitation Forcing Totals for 400-day Period 


qualitative designation 


quantitative 400-day total 


low 


32 cm 


medium 


48 cm 


high 


96 cm 


ver\' high 


160 cm 



These four complex cases were compared with the effects of "simple" constant 
evaporation and precipitation. Here the four "simple" precipitation and evaporation 
rates were specified to give constant evaporation and precipitation with the same total 
precipitation and evaporation as those of the four complex cases. These simple precip- 
itation cases were applied to the model so that the effects of discrete precipitation events 
could be compared to those of constant E-P. The resulting MED, temperature, and 
salinity time series were difierenced for comparison. The output mixed layer values of 
.MED, mixed layer temperature, and mixed layer salinity from complex and simple forc- 
ing were compared to the CTD observations. 

Advection was maintained in the model at 0.0 ppm; day. Considering the placement 
of advection in the salinity equation: 

( 3 .,, 

dt h h 

It should be noted that an increase in advection w^ould provide the same tendency for 
mixed layer salinity as would decreasing precipitation. 
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Table 3. MONTHLY PRECIPITATION AT OCEAN STATION "P" FOR 400 
DAY PERIOD 



Monthly Precipitation for 400-day Period 


month/year 


precipitation 


Dec/ 1 968 


1. 5cm 


Jan/ 1969 


0.6cm 


Feb/ 1969 


2.9cm 


Mar/1969 


3.3cm 


Apr/ 1969 


4.2cm 


May/ 1969 


3.1cm 


Jun 1969 


2.1cm 


Jul/1969 


1.3cm 


Aug 1969 


1.3cm 


Sep' 1969 


3.1cm 


Oct/ 1969 


5.5cm 


Nov; 1969 


8.0cm 


Dec 1969 


6.3 cm 



B. RESULTS AND ANALYSIS 

1. Complex Precipitation and Evaporation Forcing 

The response of the OPEL model depended upon the net precipitation for the 
cases of complex precipitation. The 400-day total precipitation was increased from low 
to ver\' high as indicated in Table 2 on page 44. An increase in total precipitation caused 
a general decrease in MLD's. as a result of increased buoyancy and decreased 
entrainment (Fig. 17 on page 47). As the precipitation level for complex precipitation 
forcing was increased from low to very high, it was larger in magnitude in winter than 
in summer, though the relative increase in summer was proportional to that of winter 
(Fig. 16 on page 46). This was due to the small amount of precipitation in summer 
compared to winter. This increase in precipitation in winter with increasing total pre- 
cipitation increased buoyancy more in winter than in summer, yielding shallow'er MLD's 
for larger complex precipitation amounts. Additionally, the largest magnitude of dilTer- 
ential MLD's occurred in winter (Fig. 17 on page 47 and Fig. 18 on page 48). 

Associated with the differential MLD flux that accompanied increased complex 
precipitation levels was a change in mixed layer temperature and salinity. Differential 
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Fig. 16. Simple and Complex Precipitation Forcing: (a) low level precipitation, 
(b) medium level precipitation, (c) high level precipitation, and (d) very 
high precipitation. 
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Fig. 17. Complex Precipitation Forcing Predicted MLD's: (a) solid line represents 
ver>’ high precipitation forced case and (b) dotted line represents low pre- 
cipitation forced case. 
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Low precipitation forcing, (a) difl'crcntial MLD, (b) difTcrcntial temper- 
ature, and (c) differential salinity^ 
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Fig. 19. Complex Precipitation Forcing Predicted Temperatures: (a) dotted line 

was low precipitation forced, (b) dashed line was medium precipitation 
forced, (c) chained-dashed was high precipitation forced, and (d) solid line 
was ver>’ high precipitation forced. 
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mixed layer temperature decreased only slightly in winter when the MLD was deep. 
Negative difTerential temperatures were induced by the winter heat loss of the shallower 
mixed layer. As the summer approached, and MLD shallowed due to increased and 
decreased winds, the differential MLD, had a more significant effect on differential 
temperature. Thus, the decreased MLD in the shallowing season concentrated more of 
the incident solar energy in a shallower layer than it would have without increased pre- 
cipitation. This produced a large positive differential temperature. 

Examination of the differential salinity showed a continual decrease in salinity 
as ver>' high complex precipitation forcing inputed more fresh water into the mixed layer. 
Deepening events which yielded differential MLD's produced differing entrainment 
events. Deeper MLD's, characteristic of lower complex precipitation forced cases, en- 
trained water of higher salinity than in the higher complex case. Larger precipitation 
events input freshwater into the surface layer inhibiting entrainment, yielding a shallower 
MLD. This reduction in MLD further concentrated the fresh water, reducing surface 
S (Fig. 20 on page 51). 

From CTD observations, actual values of MLD, temperature, and salinity were 
compared to those values predicted by the OPBL model with complex precipitation 
forcing. Verification was hindered by data voids, especially during the first 95 days of 
the study. Nevertheless, it may be concluded that the MLD's of the CTD data were best 
appro.ximated by the model using very high complex precipitation. The envelopes of the 
fluctuating MLD of all four complex precipitation forcing cases appeared to be deeper 
than those of the CTD casts. With the small number of CTD's over any synoptic pe- 
riod, an hour-by-hour comparison was not possible, although MLD was predicted to 
fluctuate diurnally by as much as 50 m'day. As the level of precipitation was increased, 
the resulting MLD became shallower. Thus, the case with the largest precipitation 
amount appeared to most closely approximate the CTD envelope for MLD (Fig. 21 on 
page 52). 

The mixed layer temperatures compared favorably with observations in all four 
complex cases. As the study period progressed, the temperatures of the actual data in- 
creased more than those from complex precipitation forcing commensurate with the ac- 
tual MLD being shallower. Thus, the case with the highest complex precipitation 
produced the greatest shallowing, yielding the highest summer temperature. The peak 
temperature from the ver\' high complex precipitation case, however, only slightly ex- 
ceeded that of the low precipitation case (Fig. 22 on page 53). 
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Fig. 20. Complex Precipitation Forcing Predicted Salinities: (a) dotted line was 
low precipitation forced, (b) dashed line was medium precipitation forced, 
(c) chaincd-dashed was high precipitation forced, and (d) solid line was 
very high precipitation forced. 
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Fig. 21. Complex Precipitation Forcing Predicted MLD and CTD MLD: (a) the 
solid line was ver\' high precipitation (breed and (b) the dotted line was 
low precipitation forced. 
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Fig. 22. Complex Precipitation Forcing Predicted and CTD Temperature: (a) 
dotted line was low precipitation forced, (b) dashed line was medium pre- 
cipitation forced, (c) chained-dashed was high precipitation forced, and (d) 
solid line was very high precipitation forced. 
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Fig. 23. Coniple.x Precipitation Forcing Predicted and CTD Salinity; (a) dotted 
line was low precipitation forced, (b) dashed line was medium precipi- 
tation forced, (c) chained-dashed was high precipitation forced, and (d) 
solid line was very high precipitation forced. 
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The model-predicted salinity values yielded less certain results. Actual CTD 
salinity data behaved in a cyclic manner during the study period. Here the initial salinity 
of 32.5 ppm immediately increased with increasing MLD. As the observed MLD shal- 
lowed, the salinity decreased, until the minimum was achieved in summer. This was 
possibly associated with shallowing due to decreased wind and precipitation. As MLD 
deepened again around day 330, salinity increased (Fig. 23 on page 54). 

Overall, this observed cyclic behavior was best approximated by modeling with 
high complex precipitation forcing. The two cases with precipitation less than the high 
case yielded salinity values that consistently exceeded those of the actual data, and were 
non-cyclic. During winter deepening for these cases, the salinity increased dramatically 
from day 0 to day 100. During spring shallowing and into summer, the salinity contin- 
ued to increase gradually, and did not exhibited a cyclical decrease. For the very high 
complex case, the salinity exhibits a decreasing trend almost immediately. Here the 
buoyant damping of the turbulence by the strong precipitation generally countered the 
shear production of turbulence attributable to the wind stress. 

Though the high complex precipitation case did agree with the cyclic behavior 
of the CTD salinity data, there were salinity observations of significant disagreement 
with the model prediction. As stated previously, peak salinity in winter w'as less than the 
actual value. In summer the salinity minimum w'as larger than the observed minimum 
by ~ 0.7 ppm. 

2. Constant Precipitation and Evaporation Forcing 

The results obtained during this study were similar to those for the complex 
precipitation forcing study. Upon comparing the various magnitudes of simple precipi- 
tation forcing, it w'as expected that increasing levels of constant E-P would have a 
magnifying effect on the mixed layer parameters. For simplicity, high and low^ precipi- 
tation forcing cases w'ere compared here. The increased level of precipitation in the high 
case yielded a shallow'er overall MLD than the low^ case (Fig. 24 on page 56). This 
shallowing effect was stronger in winter than in summer, as MLD's tended to be deeper 
in winter. The peak predicted in the shallower period in the low case (Fig. 24 on page 
56 point a) did not appear in the high case. 

It appeared that winter temperatures were w^ell verified for all simple and com- 
plex forced cases by the actual mixed layer temperatures, as the increased precipitation 
yielded only slightly more shallowing. In summer, how'ever, differential temperatures 
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Fig. 24. Simple Precipitation Forcing Predicted MLD'S: (a) the solid line was 
very high precipitation forced and (b) the dotted line was low precipitation 
forced. 
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Fig. 25. Simple Precipit.itioii Forcing Predicted Temperature: (a) dotted line was 
low precipitation forced, (b) dashed line was medium precipitation forced, 
(c) chained-dashed was high precipitation forced, and (d) solid line was 
ver>’ high precipitation forced. 
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increased as the high case absorbed more solar heat in a shallower layer than the low 
case. 

As with the complex precipitation cases, the constant low and medium precipi- 
tation predictions yielded a proportional decrease in salinity. The high precipitation case 
proved to be the most cyclic in time over the year. The very' high case yielded an almost 
monotonic decrease in salinity as a result of the large constant surface fresh water influx 
(Fig. 26 on page 59). 

From comparing simple to complex precipitation forcing, it was noted that the 
complex case produced relatively more precipitation in winter than in summer, and there 
was greater evaporation in summer than in winter. The constant E-P forcing effect of 
the simple case thus had a larger precipitation rate during periods of no rain for the 
complex case (in winter when low winds occur and during the summer). In winter dif- 
ferential .VIED fluctuated both positively and negatively about zero. It appeared to 
fluctuate negatively during a strong complex precipitation event yielding shallowing. It 
fluctuated positively during non-rain periods for the complex case, producing relative 
shallowing in the simple case. In winter these differential MLD values exceeded those 
of summer. 

Differential temperature in winter showed approximately equal temperatures for 
the four rain levels in both the simple and complex cases. Additionally, any fluctuation 
in MLD had a relatively small effect on the differential temperature in the winter season. 
In summer a negative differential temperature was induced by precipitation being greater 
in the simple case than in the complex case. 

With respect to salinity, differential salinity tended to be positive in winter. 
Here, the salinity of the complex case was usually less than that of the simple case due 
to the greater amount of precipitation in the complex case. It did fluctuate during this 
period, as periods of low winds yielded a decreased complex precipitation allowing sim- 
ple precipitation to exceed it and produce a negative fluctuation in differential salinity. 
In summer a positive d’flerential salinity pre-dominated. Here, the lower precipitation 
for the complex case allow'ed the salinity of this case to increase. 

The very’ high simple precipitation case provided the best approximation to ac- 
tual .MLD of all study cases, both simple and complex. The large amount of precipi- 
tation produced consistently shallower MLD's than the other simple cases. Additionally, 
it provided shallower .MLD's than occurred during the complex case when forced by low' 
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Fig. 26. Simple Precipitation Forcing Predicted Salinity: (a) dotted line was low 
precipitation forced, (b) dashed line was medium precipitation (c) 
chained-dashed was high precipitation (d) solid line was very high precip- 
itation. 
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forcing output by Low precipitation forcing, (a) differential MLD, (b) 
differential temperature, (c) differential salinity, and (d) complex precipi- 
tation forcing. 
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Fig. 28 . Simple, Complex, and Actual Mixed Layer Values: In all three figures 

the solid line represents simple high forcing, the dotted line represents 
complex forcing. 
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winter winds and in all summer cases when low winds occur and MLD's tend to be 
shallowest (Fig. 28 on page 61). 

In winter, there was close agreement between all complex and simple forced 
cases and actual temperature. In summer, the shallowest MLD's of the simple forced 
case produced the highest temperature, which though best approximated the actual 
temperature in summer, were still less than that of the actual case by 2.4°C. 

As with complex forcing, the high simple case provided the best approximation 
to the actual salinity values. The simple high rain case yielded a larger negative differ- 
ential salinity in winter and a low’er salinity in summer than the high complex case, while 
still retaining a cyclic behavior with a one-year period. Its cyclic response was not as 
clear as in the complex case, however. 
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VI. SUMMARY 



This study was conducted to examine the discrete precipitation event effects on the 
short-term and seasonal evolution of ocean mixed layer temperature and salinity struc- 
ture. Data for this study was for the years 1968 and 1969 at Ocean Station "P" located 
in the Northeast Pacific Ocean. Two kinds of experiments were performed; 

1. the first was to simulate the response of the mixed layer to discrete (isolated) pre- 
cipitation events, and 

2. the second was a quasi-realistic "complex" precipitation forcing experiment, re- 
presenting a realisnc distribution in time of synoptic events over a 13-month pe- 
riod. 

The hypothetical experiment was conducted using twelve separate single discrete 
precipitation events over a 400-day period. It examined the specific effects of a discrete 
precipitation event as a function of season without contamination by other complex 
rainfall history. This experiment showed that the effect of a single discrete rain event 
can vaiy with season and the MLD (at onset of the rain event), as well as with the in- 
tensity of the wind stress. Additionally, it showed that a single rain event can have a 
lasting effect on mixed layer depth and temperature for up to 55 days after the event 
stops. 

For the main experiment a "complex" precipitation modeling scheme, producing 
discrete daily precipitation events from surface meteorological observations was devised 
and tuned to yield realistic precipitation input for the mixed layer model. The precipi- 
tation model was adjusted to give differing amounts of total precipitation. Four different 
precipitation "intensities" were tested: 32 cm'400-day, 48 cm; 400-day, 96 cm, 400-day, 
and 160 cm, '400-day. The 50 cm'400-day case was chosen to yield the same net precip- 
itation as the observed precipitation data. These quasi-realistic precipitation forcings 
were applied to the mixed layer model. The resulting values of MLD, temperature, and 
salinity for this 400-day period were examined by comparing them with output values 
when the model was forced by "simple" constant net precipitation with the same total 
precipitation as in the complex cases. Finally, the values predicted by the model under 
various precipitation forcing conditions were compared to actual CTD values. 

In all study cases the MLD's approximated the MLD cycle of the actual data. The 
observed MLD's were shallower than predicted by the model by 10m or more. Thus, the 
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model runs with the most precipitation best approximated the observed MLD's. It was 
noted that the greatest discrepancy occurred in the first winter season of the 400-day 
study period. However, a large data void existed during this time period. 

As with MLD, the temperatures of these studies generally approximated the tem- 
perature cycle and form of the actual data. The temperature in summer provided the 
greatest discrepancy in temperature magnitude with actual temperatures (exceeding 
those of the output data by at least ~ 2.5° C in the complex case and ~ 2.5°C in the 
simple case). 

Predicted salinity provided the greatest discrepancy between the modeled and the 
observed cycle. Rain rates of 32 cm/yr and 48 cm'yr (that equal to the observed total) 
yielded salinity values which were non-cyclic and, overall, higher in value than the ob- 
servations. The 96 cm '400-day case (both complex and simple) best approximated the 
salinity observations. Predicted salinity was less than observed in winter (by ~ O.lOppm 
in the complex case and ~0.11ppm in the simple case) and greater than observed in 
summer (by ~ 0.01 ppm in the complex case and ~ 0.01 ppm in the simple case). The 160 
cm 400-day precipitation cases produced a constant freshening of the mixed layer in 
both winter and summer. 

In the final analysis, while the discrete precipitation events had a different effect on 
the mixed layer than the simple cases, the difference was slight. Increased levels of pre- 
cipitation did enhance the difference between constant rain and the quasi-realistic rain 
cases, as seasonal differences in precipitation were amplified. The 96 cm/400-day 
quasi-realistic precipitation forcing case was selected as the most likely precipitation 
forcing for Ocean Station "P" during this time period. Despite magnitude discrepancies 
with observed rain at Station "P," the form and cycle of the MLD's and temperatures 
with respect to time were well approximated by this level of precipitation. The most 
important point for selection of this precipitation case was, however, the agreement with 
the actual salinity cycle. Here, the simple high precipitation case also approximated the 
shape of the salinity curve, but was not as cyclical, especially during the winter in 1969. 
This was probably because the complex precipitation case produced some periods with- 
out rain allowing for shallowing, and, therefore, decreased salinity. Even in t^inter, while 
the simple cases produced continuous rain. One such realistically modeled period acted 
to enhance cyclical behavior of salinity in the early winter of 1969, adding to MLD 
deepening and adding less fresh water during this period than that of the simple case. 
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There are several reasons why observed mixed layer temperature and salinity values 
may not agree with model predicted values: 

1. the inaccuracies and voids in CTD data, 

2. the inaccuracies in model forcing, especially in producing realistic complex precip- 
itation forcing, 

3. the inaccuracies in model initial conditions for T(z) and S(z), 

4. the inaccuracies in representing advection, and 

5. the mixed layer theory inaccuracies. 

Data voids do exist in the available CTD data at Ocean Station "P." The largest 
data void occurred during the first 95 days of the study. At this time, there w^as an ob- 
servation at day 0 (Xov 23, 1968) and then no observations until the 95th day (Feb 26, 
1968). Additionally, data w'as sparse from days 95 to 180 and days 300 to 400. This 
would not be significant for data w'ith small variance, as with the sea surface temper- 
ature. How'ever, for data with strong variation such as the .MLD and salinity data, the 
resulting analysis may be biased by lack of information at times of critical changes in 
mixed layer evolution. 

Inaccuracies in prec pitation forcing for the model may exist. Four net precipitation 
rates, each multiples of the first, were chosen to bracket the most realistic precipitation 
amounts for Station "P." The fact that the 48 cm/400-day case did not approximate the 
salinity cycle as well as the 96 cm/400-day case, suggests that the observed precipitation 
rate may be low'er than the actual amount. Allen (1963) outlined the significant number 
of potential rain error measurement sources onboard ship at Ocean Station "P." These 
were: 

1. the wind effects, 

2. the shadowing by nearby objects, 

3. the evaporation in the gauge, 

4. the sea spray entering the gauge, and 

5. the tilting of the gauge orifice due to roll and pitch of the ship. 

Additionally, the rain model scheme used in this study was very crude. It may have 
inaccurately approximated the strength and duration of storm events, and would have 
been completely insensitive to local convection events (rain w'ithout wind). Such local 
convective events may account for a lower observed salinity, shallower MLD's, and 
higher temperatures in summer. 
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Inaccuracies in meteorological data used in this study to compute heat and mo- 
mentum fluxes may also exist. This data was collected in 1968 and 1969, and collection 
techniques have improved since. In addition to the wind stress and heat flux errors in- 
troduced directly into the model by this data, these errors could produce unrealistic 
precipitation and evaporation events as well. 

Predicting horizontal advection is diflicult, and given only surface meteorological 
data and scattered CTD data it is not possible with any accuracy. Advective events have 
been known to occur at this station, as discussed in Chapter I. Such an event could 
greatly modify mixed layer structure. Advection was neglected for these studies, recog- 
nizing that an increase in advection of water of higher salinity has the same effect as 
decreasing rain and is necessary to balance the salinity budget in the longer interannual 
term. This was another reason for examining various rain amounts. 

It is recommended that this experiment be repeated using observed discrete precipi- 
tation if possible. Alternatively, Tucker (1961) suggested a rainfall prediction scheme 
utilizing weather codes, using the standard international two-digit system. These 
follow-on studies could help answer the question as to the discrepancy in actual versus 
measured rainfall. Finally, precipitation from such sources as satellite imaging could be 
interpreted as with Hakkarien and Adler (1987) and applied to the mixed layer model. 
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